########################
###Summary Statistics###
########################
library(foreign)

# set working folder. 
setwd("~/OneDrive - Indiana University/FromGoogle/syd/")
# read dataset
full.2003.2018 <- read.csv("modmeddataset7225.csv")

#Cancel sci notation
options(scipen = 999)
###clear_zoon_confirm
c(min(full.2003.2018$clear_zoon_confirm, na.rm = TRUE),median(full.2003.2018$clear_zoon_confirm, na.rm = TRUE),mean(full.2003.2018$clear_zoon_confirm, na.rm = TRUE),max(full.2003.2018$clear_zoon_confirm, na.rm = TRUE),sd(full.2003.2018$clear_zoon_confirm, na.rm = TRUE))
###ag_crops
c(min(full.2003.2018$ag_crops, na.rm = TRUE),median(full.2003.2018$ag_crops, na.rm = TRUE),mean(full.2003.2018$ag_crops, na.rm = TRUE),max(full.2003.2018$ag_crops, na.rm = TRUE),sd(full.2003.2018$ag_crops, na.rm = TRUE))
###Forest.coverage_bin
c(min(full.2003.2018$Forest.coverage_bin, na.rm = TRUE),median(full.2003.2018$Forest.coverage_bin, na.rm = TRUE),mean(full.2003.2018$Forest.coverage_bin, na.rm = TRUE),max(full.2003.2018$Forest.coverage_bin, na.rm = TRUE),sd(full.2003.2018$Forest.coverage_bin, na.rm = TRUE))
###NDVI_mean_anom
c(min(full.2003.2018$NDVI_mean_anom, na.rm = TRUE),median(full.2003.2018$NDVI_mean_anom, na.rm = TRUE),mean(full.2003.2018$NDVI_mean_anom, na.rm = TRUE),max(full.2003.2018$NDVI_mean_anom, na.rm = TRUE),sd(full.2003.2018$NDVI_mean_anom, na.rm = TRUE))
###NDVI_mean_anom
c(min(full.2003.2018$fires.cell.count_anom, na.rm = TRUE),median(full.2003.2018$fires.cell.count_anom, na.rm = TRUE),mean(full.2003.2018$fires.cell.count_anom, na.rm = TRUE),max(full.2003.2018$fires.cell.count_anom, na.rm = TRUE),sd(full.2003.2018$fires.cell.count_anom, na.rm = TRUE))
###ag_crop
c(min(full.2003.2018$ag_crop, na.rm = TRUE),median(full.2003.2018$ag_crop, na.rm = TRUE),mean(full.2003.2018$ag_crop, na.rm = TRUE),max(full.2003.2018$ag_crop, na.rm = TRUE),sd(full.2003.2018$ag_crop, na.rm = TRUE))
###Forest.coverage_bin
c(min(full.2003.2018$Forest.coverage_bin, na.rm = TRUE),median(full.2003.2018$Forest.coverage_bin, na.rm = TRUE),mean(full.2003.2018$Forest.coverage_bin, na.rm = TRUE),max(full.2003.2018$Forest.coverage_bin, na.rm = TRUE),sd(full.2003.2018$Forest.coverage_bin, na.rm = TRUE))
###log.nl
c(min(full.2003.2018$log.nl, na.rm = TRUE),median(full.2003.2018$log.nl, na.rm = TRUE),mean(full.2003.2018$log.nl, na.rm = TRUE),max(full.2003.2018$log.nl, na.rm = TRUE),sd(full.2003.2018$log.nl, na.rm = TRUE))
###log pop
c(min(full.2003.2018$logged.pop, na.rm = TRUE),median(full.2003.2018$logged.pop, na.rm = TRUE),mean(full.2003.2018$logged.pop, na.rm = TRUE),max(full.2003.2018$logged.pop, na.rm = TRUE),sd(full.2003.2018$logged.pop, na.rm = TRUE))
###spei
c(min(full.2003.2018$spei, na.rm = TRUE),median(full.2003.2018$spei, na.rm = TRUE),mean(full.2003.2018$spei, na.rm = TRUE),max(full.2003.2018$spei, na.rm = TRUE),sd(full.2003.2018$spei, na.rm = TRUE))
###p_anom
c(min(full.2003.2018$p_anom, na.rm = TRUE),median(full.2003.2018$p_anom, na.rm = TRUE),mean(full.2003.2018$p_anom, na.rm = TRUE),max(full.2003.2018$p_anom, na.rm = TRUE),sd(full.2003.2018$p_anom, na.rm = TRUE))
###t_anom
c(min(full.2003.2018$t_anom, na.rm = TRUE),median(full.2003.2018$t_anom, na.rm = TRUE),mean(full.2003.2018$t_anom, na.rm = TRUE),max(full.2003.2018$t_anom, na.rm = TRUE),sd(full.2003.2018$t_anom, na.rm = TRUE))
###spllag_zoon
c(min(full.2003.2018$spllag_zoon, na.rm = TRUE),median(full.2003.2018$spllag_zoon, na.rm = TRUE),mean(full.2003.2018$spllag_zoon, na.rm = TRUE),max(full.2003.2018$spllag_zoon, na.rm = TRUE),sd(full.2003.2018$spllag_zoon, na.rm = TRUE))
###spllag_ndvi
c(min(full.2003.2018$spllag_ndvi, na.rm = TRUE),median(full.2003.2018$spllag_ndvi, na.rm = TRUE),mean(full.2003.2018$spllag_ndvi, na.rm = TRUE),max(full.2003.2018$spllag_ndvi, na.rm = TRUE),sd(full.2003.2018$spllag_ndvi, na.rm = TRUE))
###spllag_fires
c(min(full.2003.2018$spllag_fires, na.rm = TRUE),median(full.2003.2018$spllag_fires, na.rm = TRUE),mean(full.2003.2018$spllag_fires, na.rm = TRUE),max(full.2003.2018$spllag_fires, na.rm = TRUE),sd(full.2003.2018$spllag_fires, na.rm = TRUE))




